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Transverse Stress Decay in a Specially Orthotropic 
Strip under Localized Normal Edge Loading 

W. B. Fichter 
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ABSTRACT 

Solutions are presented for the stresses in a specially orthotropic infinite strip which 
is subjected to localized uniform normal loading on one edge while the other edge is 
either restrained against normal displacement only, or completely fixed. The solutions 
are used to investigate the diffusion of load into the strip and in particular the decay 

of normal stress across the width of the strip. For orthotropic strips representative of a * 

broad range of balanced and symmetric angle-ply composite laminates, minimum strip 

widths are found that ensure at least 90% decay of the normal stress across the stop. In 

addition, in a few cases where, on the fixed edge the peak shear stress exceeds the normal 

stress in magnitude, minimum strip widths that ensure 90% decay of both stresses are 

found. 

To help in putting these results into perspective, and to illustrate the influence ot 
material properties on load diffusion in orthotropic materials, closed-form solutions for 
the stresses in similarly loaded orthotropic half-planes are obtained. These solutions 

are used to generate illustrative stress contour plots for several representative laminates. i 

Among the laminates, those composed of intermediate-angle plies, i.e., from about 30“ 
to 60°, exhibit marked changes in normal stress contour shape with stress level. 

The stress contours are also used to find 90% decay distances in the half-planes. 

In all cases, the minimum strip widths for 90% decay of the normal stress exceed the 
90% decay distances in the corresponding half-planes, in amounts ranging from only 
a few percent to about 50% of the half-plane decay distances. The 90% decay dis- 
tances depend on both material properties and the boundary conditions on the supported 
edge. 

The widespread and growing use of composite materials in recent years has prompted 
increased interest in analysis of the stresses and deformations in elastic anisotropic materials. 

As a result, it is well known that stress decay rates in anisotropic elastic bodies can differ 
markedly from those in similarly loaded isotropic bodies. (See, for example, [1] and [2], 
where numerous studies are cited.) Many of the cited studies have focused on questions 
concerning the applicability of St. Venant’s pnnciple to anisotropic bodies, particularly 
elastic strips that are loaded in the longitudinal direction. One finding of such investigations 
is that the distance over which end effects in axially loaded anisotropic strips diffuse and 
decay to insignificant levels can differ by an order of magnitude or more from that in 
an isotropic strip, a result which has obvious implications in the determination of elastic 

properties of composite materials from laboratory tests of small specimens. In [3], solutions , 

were found for the stresses in a specially orthotropic half-plane subjected to localized , j 

i ■ * 
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£,, £ v , G, v 

H i v ■* Hv.c 

h 


k , /c i , ^2 » q< r 
p 


Nomenclature 

coefficients in solution to m, u 
transformed governing 
equation 

elastic constants for 
orthotropic materials 
dimensionless width 
of strip 

minimum strip width 
required to ensure 
90 % decay of o\ 
characteristic parameters 
magnitude of applied 
normal traction 


x , v 

y\i 

y 

£ x * £ y * T vv 


cr t , cr v , T t 
<P(JC. y) 


in-plane displacements 
rectangular Cartesian 
coordinates 

90 % decay distance for 
0 V in the half-plane 
characteristic root of 
transformed 
governing equation 
elastic strains 
Fourier transformation 
variable 
elastic stresses 
Airy stress function 


Subscripts C and 5 indicate Fourier sine and cosine transforms, respectively. 


self-equilibrated edge loading, and a simple approximate formula was found for a lower 
bound on ihe normalized 90% decay distance for the normal stress. 

There are practical concerns as well for the decay of stresses across the width of elas 
strips. Many modem aircraft and spacecraft stractnral components 

edee are no more than a specified fraction of the applied loads. 

8 Making this tvpe of information available to designers is the primary poipose of thl 
article A specially orthotropic (principal axes of orthorropy aligned with coordinate axes) 
dast smpTsut ec.ed to legalized uniform normal loading on one edge, while the opposne 
edge is drtier restrained against normal displacement only, or complerely fixed. For bo,h 

similarly loaded specially orthotropic and isotropic half-planes. 

ANALYSIS 

The specially orthotropic strip studied here is defined by -oo < x< oo. 0 _^> - ^ 
The engineering constants are E x , E y . G x> , Tty an ri vt . w ere < >* 
stress conditions the stresses, strains, and displacements are related by 

e* = u.x = t 1 - J- £v = r, v = £7 ~ E < 


E x ■” E v 

Ttv 

Ytv = Tl.v + l ’.Jt = 

V'JtV 


( 1 ) 


The stresses are 


defined in terms of the Airy stress function by 


ff, = <P.v 


CTv = <P.xx 


Ttv = -<P..c 


( 2 ) 
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where the comma followed by a subscript indicates partial differentiation with respect to the 
subscript variable. Thus, the plane-stress equilibrium equations are satisfied automatica y, 
and the compatibility equation can be written as 


YV.xxxx + - 2n,v)<P «yy + (p w- v " ° 


(3) 


Two boundary-value problems for the specially orthotropic strip, referred to herein 
as p^blems A^d B, ar/analyzed here. In both problems, uniform ^nonma. .ac on 
■ i- j f rPP pHge v = 0 over an interval of lengtn z, wmcn, iui 

simplicity is centered on the v axis. (Note that the stresses in all the solutions presented 
here can^be easily converted to solutions for applied traction over an interval of lengt 

f-dV are «. e<,ual » zero. The sohmon ro problem A . presenred 

first. 


Problem A— Partially fixed upper edge 
The boundary conditions are 


T,. v (.t, 0) = 0 


\p |je| < 1 

"' <j:0)= |o W>1 

T tv (.r, h) = 0 
v(x, h) = 0 


The unbounded domain and the discontinuous edge loading suggest a 
approach. The Fourier cosine and sine transforms are defined by 


<Pc(3^ y) 
(p 5 (A, y) 



<p(x , y) cos \x dx 
cp(x, y) sin Xx dx 


or 


and inversely. 


cp(x, y) 
cp(x, v) 



cp c (A, y) cos Xx dX 


i p s (A, y ) sin \x dX 


or 


(4) 

(5) 

( 6 ) 
(7) 


Fourier transform 


( 8 ) 


( 9 ) 
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In terms of the transformed stress function, the transformed stresses and displacements 
are 


®xC — <Pc,yy tfyC = “ ^~<Pc Tty 5 = ^<PC\y 

“s = (^-J<Pc.w + and (10) 

yc = (a I £;) <PC vvv + ( £7 “ G^) <PC v 

which anticipate obvious symmetries about the v axis in both problems. Application of the 
appropriate transforms to the governing equations yields 


<PC.yyyy ~ ~ 2q, v ^ A 2 <p C .yy + ^|^A 4 <PC = 0 


along with 


cpc.y(A.O) = 0 


, p2 ( sin A A 

<Pc(A, 0) = - Py j- (^— \ = -/(A) 


<Pc.v( A. h) = 0 


and 


<Pc.y»(A, h) - {^r~ - Vry^A 2 <p c .y(A, h) = 0 

which in view of Eq. (14) simplifies to 


( 11 ) 


( 12 ) 

(13) 

(14) 


<PC,y>y(A. h) = 0 (15) 

When the solution to Eq. ( 1 1) is assumed to have the form <p c = F( A)e* Yy , the following 
characteristic equation is obtained: 

t4 -(^- 2,, ") t!+ (|) =o ,i6 ' 

With the definitions 2 r = (E x /G xy — 2r| xy ), k 4 = E x / E yy Eq. (16) becomes 

Y 4 -2ry 2 +* 4 = 0 (17) 

The form of the solution for <pc(A, y) depends on the nature of the characteristic 
roots, which depends on the material constants through the parameters r and k , as seen 
in Eq. (17). With the physically meaningless special cases disregarded, the three material 
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types of engineering interest are: 

Type /: 0 < k 2 < r. Let q = /P- - EJE y , and k { = y/FTq, k 2 = y/F^q. Then 
Y = ±& i, ±k 2 (real, distinct). 

Type //: 0 < k 2 = r. Then y = ±k , ±k (real, r epeated). 

Type 111 : 0 < \r\ < k 2 . Let k { = y/Uc 2 + r)/2, k 2 = y/ik 2 - r)/2. Then y = ±(*, ± ik 2 ) 
(equal and opposite conjugate pairs). 

Most balanced and symmetric angle-ply laminates are of type I or type III. However, one 
type II material of great practical interest is the isotropic (or quasi-isotropic) material. The 
isotropic solution can be found from a type I or type III solution through an appropriate 
limiting process, as well as by dealing directly with the isotropic formulation. 

Solution for type / material 

The transformed stress function has the form 

<p c (A, y) = A(A)cosh A£i_y + fl(A)sinh A*,y + C(A) cosh A/by + D(A)sinhM 2 y (18) 

Substitution of Eq. (18) into the transformed boundary conditions, Eqs. (12M15), leads to 
simultaneous algebraic equations to be solved for A, B, C, and D, which are then substituted 
into Eq. (18) to completely define cp c (A. y ). Then the use of Eqs. (8H 10) gives for the 
stresses at any point in the strip, 


poo 

cr t = (2pjn)k\ki I [£j sinh Xk 2 h cosh Xk\(h — y) 

Jo 

- ko sinh Xk\h cosh A k 2 (h - y)] sin A cos Xx d\/A 

poo 

cr v = ( 2 p/tt) I [k] sinh \k[h cosh \k 2 (h - y) 

Jo 

- k 2 sinh A k 2 h cosh A k\{h - y )] sin A cos XxdX/A 

poo 

Try = (• 2 p/ir)k\k 2 I [sinh \kih s'mh\k 2 (h — y) 

Jo 

- sinh \k 2 h sinh \kfh — y)] sin A sin Ax d\/A 

where 

A = A(£| sinh \k[h cosh \k 2 h - k 2 cos,h\k x h sinh Xk 2 h) 
On the supported edge y = h, the nonzero stresses are 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


<J x {x,h) = 2{k\ + kl)p/7T / [^! sinhA^A - k 2 sinh A£i/i]sin Acos Ax rfA/A (23) 

Jo 

Gy(x, h) = (2p/:i) f [£t sinh Xk\h — k 2 sinh Xk 2 h] sin A cos Xx dX/XA (24) 

Jo 

Solution for type III material 

With the characteristic-value definitions y = k\ + ik 2 and y = k\ — ik 2 for type III 
materials, the transformed stress function has the form 


cp c (A, v ) = A( A) cosh Ayv + B( X) sinh Ayv + C( A) cosh Ay y + D( A) sinh Ay v (25) 


158 


W. B. Fichter 


Satisfaction of the transformed boundary conditions, Eqs. (12H15), and the subsequent 
use of Eqs. (25) and (8)-(10), yields the following expressions for the stresses at any point 

in fnp ctrtn* *' * 


CT. t — 2p(k\ +kl)/i t I [F| coshM,(/i - y)cosXk 2 (h - y) 

J 0 

+ F 2 sinh Xk;(h - y) sin Xk 2 (h - y)] sin A cos Ax dX/A 

where 


F\ — k 2 sinh A k\h cos A k 2 h — k[ cosh Xk\h sin Xk 2 h 
F 2 = k 2 cosh Xkih sin Xk 2 h + k { sinh Xk t h cos Xk 2 h 


/*O0 

Oy=2p/n I (& 2 [sinh Xkih cos Xk->y cosh A/ti 

Jo 


y ) 


(27) 


(28) 


(29) 


(h - y) 

JO 

- sin Xk 2 h sinh A k x y sin A k 2 (h - y)] + [sin Xk 2 h cosh Ait, y cos Xk 2 (h - 
+ sinh A k, h sin A k 2 y sinh A (h - y)]J sin A cos A* dX/ A 

n oo 

T xy = ^p{k\ + k^)Jn J [sinh Xk[h sin Xk 2 y cosh Xk\(h — y) 

~ sin Xk 2 h sinh A k { y cos A k 2 (h - y )] sin A sin Ax dX/A 

where 

A = X(k 2 sinh Xk\h cosh Xk\h + k\ sin A k 2 h cos Xk 2 h) 

On the supported edge v = h, the nonzero stresses are 

h) = 2p{k] + kj)/7T f [k 2 sinh Xk\h cos Xk^h 
Jo 

~~ k\ cosh Xk\h sin A k 2 h ] sin A cos Xx dX/ A 

roc 

<J y (x, h) = 2p/n I [k 2 sinh Xk[h cos Xk 2 h 
Jo 

+ k { cosh Xk[h sin AX: 2 ^] sin A cos Ax dX/A 
Solution for isotropic material 

The details of the solution are presented in the Appendix. The formulas for the stresses 
at any point in the strip are 


(30) 


(31) 


/•00 

o x = 2p/n / [sinh Ay + sinh A(2A - y) - X{2h - y) cosh Ay 
J 0 

- Ay cosh A(2 h - y)] sin A cos Ax dX/A 

r oo 

°y = 2 P/ n / t sinh \y + sinh X(2h - y) + A(2 h - y) cosh Ay 

J 0 

+ cos h X(2h - y)] sin A cos XxdX/A 

T r 


(32) 


(33) 




r oo 

^P/ 71 sin h — y) - X(2h - y) sinh Ay] sin A sin Ax dX/A (34) 
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where 


A = A(sinh2A/i + 2A/i) 


(35) 


On the supported edge y — h, the nonzero stresses are 

o’.tCv, h) — 4 p/n s [sinh A h - A/i cosh A h] sin A cos A.c d\f A (36) 

/»oo 

OyCt, h) — 4p/n I [sinh Xh + Ah cosh A/ 2 ] sin Acos Ajc dA/A ( 37 ) 

J 0 


Problem B — Fixed upper edge 
The boundary conditions are 


T ty (x, 0 ) = 0 
v y (x. 0 ) = 
u(x, h) = 0 


P 1*1 < 1 

0 U| > 1 


and 


u(.t, h) = 0 

In terms of <p c , the transformed boundary conditions are 

<Pc.v(A, 0) = 0 


<p c «.o)=- P yi(^).- /( x, 


(pc. y .v(A, h) + A 2 r|, y <p c (A. h) = 0 


and 


A) + A 




(Pc, v (A, h) = 0 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


Solution for type / material 

The use of Eq. (18) in conjunction with Eqs. (42H45) yields a set of simultaneous 
equations which is solved for A, B , C, and D. These are used along with Eqs. (18) and 
( 8 H 10 ) to obtain the expressions for the stresses at any point in the strip: 

/* OO 

<J x =2pk\k 2 /n {(£,** - r\ 2 )(k^ cosh M { y + k\ cosh A£ 2 .v) 

+ (£ + h) [^ 1^2 sinh \k 2 h sinh AA: i (/i — v) — k\ cosh \k\h cosh \k 2 (h — v)j 

-h (i - n) 2 [*i *2 sinh \k { h sinh A k 2 (h - y) - k\ cosh \k 2 h cosh A k x (h - y)]) 
sin Acos A.r d\/A 


( 46 ) 
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(j v = 2p/7z I — f|) 2 cosh A/c 2 ^ cosh A/cj (/i — y) 

Jo 

+ (£. + r|) 2 cosh Xk\h cosh A k 2 (h — y) — (e — r|)(cosh Xk { y + cosh A/^y)] 

— k\(e — T ]) 2 sinh Xk[h sinh A^(^ — y) — kj(t + r|) 2 sinh A k 2 h sinh A k\(h — v)} 
x sin A cos Ax dXj A (47) 

f°° j 

t xv — 2pk\k 2 jn I j(£, 2 — r\ 2 )(k\ sinh Xk[ v 4- k 2 sinh Xk 2 y) 

Jo 

+ k i(£, - r|) 2 [cosh A^/i sinh Xk[(h — >’) — sinh Xk\h cosh Xk 2 (h — y)] 

+ kz(L + rj) 2 [^osh A k\h sinh A k 2 (h — v) — sinh A k 2 h cosh A k[(h — y ) ] } 

sin A sin Xx dX/A (48) 

where 


l = 1/2 (Jtf + kj + 2v X y) n = l/2(* 2 - k\) (49) 

A = A{2ifcifc 2 [(f. 2 +n 2 ) cosh AA:,/i cosh Ai : /i - (£, 2 - r| 2 )] 

— [fc 2 (E, — r|) 2 + k\{l, + r|) 2 ] sinh Xk\ h sinh Xkih } (50) 

Along the fixed edge y — h y the formulas simplify to 

d x (x, h) — r\ xy <y y (Xy h) (as required by the fixed-edge condition) (51) 

roc 

c 7 y (x, h ) = 4pk\k 2 r\/7X I [(E, + rj)cosh Xk\h — (£ — r|)cosh Xk 2 h] sin A cos Xx dX/A (52) 
Jo 

fOC 

Tj Ky (x, h) = 4pk\k 2 V[/n J [k\ (E, - r\) sinh Xk [h—k 2 (^ + r\) sinh Xk 2 h] 

Jo 

sin AsinAxdA/A (53) 


Solution for type III material 

Substitution of Eq. (25) into the transformed boundary conditions (42H45) produces 
a set of simultaneous equations which are solved for A, B , C\ and D. Substitution of their 
values into Eq. (25) and subsequent use of Eqs. (8) — (10) yields for the stresses, 


a x = 2yy p/tt I Re[ggy 2 cosh Ayy + yyg 2 sinh Xyh sinh Ay (h — y) 

Jo 

— y 2 g 2 cosh Xyh cosh Ay (h — y )] sin A cos Xx dX/ A (54) 

<7 y = 2 p/n I Re[yyg 2 cosh Xyh cosh Ay (h — v) — gg cosh Ayy 

jo 

— y 2 g 2 sinh Xyh sinh Xy(h — y)] sin A cos Ax dX/A (55) 

? 

t xy = 2yyp/iz I Re[yg cosh Ay/i sinh Ay (/i — y) + ygg sinh Ayy 
Jo 

— yg 2 sinh Xyh cosh Xy(h — y )] sin A sin Ax d A/ A (56) 


where again y = k\ 4- ik 2 , y ~ k { ~ ik 2 , and g = y 2 + t| T> „ g = y 2 4- v xy9 
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Re{c) indicates the real part of the complex number z y and 

A — A{yy[Re(g 2 )cosh Ayh cosh Ay h — gg] - Re(y 2 g 2 ) sinh Ay/t sinh Ay/t} 

or, in explicitly real terms, 


A = A((fc 2 -h k\)\Ak\k\ — (k 2 — k\ + ri Jtv )‘](sinh 2 Ak[h 4- cos 2 A&2^) 

+ (*? + *I)[(*f - ^*2 + n xy ) 2 + 4kjkl] +{(*?- k 2 )[{k 2 - k\ + Tlx ,) 2 - 4*r*2 2 ] 
4- 8 k\k\{k\ — k\+ rj JC> ,)}(sinh 2 Ak\h 4- sin 2 Akih)) (57) 

On the fixed edge y = h, the expressions for the stresses become 

<j x {x , h) = r\ xv a y (x J h) (as required by the fixed-edge condition) (58) 

/•OO 

& y (x, h) = %k[ki{k\ 4- k\)p/n I [2k \ki cosh Ak\ h cos Akih 

Jo 

4- (k 2 — k\ 4- r[ xy ) sinh Ak\h sin A/^] sin A cos A* dA/A (59) 


and 


r oc 

T ^(jc , /i) = 8 /:i^ 2(^? + fcf) /7 / 7t J [&i (A: 2 4- & 2 + !4t V ) cosh A£[/i sin A&2^ 

— + k\ — \i xy ) sinhA/^/icos A/^] sin A sin Ax dA/A (60) 


Solution for isotropic material 

The details of the solution are presented in the Appendix. The formulas for the stresses 
at any point in the strip are 

a x = 2p/n f {(3 — i>) sinh Ah [sinh A (h — y ) — Ay cosh A(h — y )] 

Jo 

+ [4r|/(l 4-n)-(l 4- r\)A 2 h(h - y)] cosh Ay 

+ [(1 4- r|)A/i 4- 2Ay] sinh Ay} sin A cos Ax dA/A (61) 

r 

<s y = 2 p/n I {(3 — rj) sinh A/z(sinh A (h — y ) 4- Ay cosh A(h — y)) 

Jo 

4-[(l 4“ r\)A 2 h(h - y) 4- 4/(1 4 - t])] cosh Ay 

4- [(1 4- r\)Ah — 2Ay] sinh Ay} sin A cos Ax dA/A (62) 

roo 

x xy —2pjn I {(3 — r|)Ay sinh Ah sinh A(h — y) 4- 2Ay cosh Ay 

Jo 

- [2(1 — rj)/(l + r|) 4- (1 4- r\)A 2 h(h - y)]} sin A sin Ajc dA/A (63) 


where 


A = M(3 - n) sinh 2 A/i + (1+ r\)X 2 h 2 +4/(1+ n)] 


(64) 
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On the fixed edge y = h. 


<r t (x, h) = r]CF v (jt, h) (as required by the fixed-edge condition) 

poo 

cr v .(.r, h) = 2p/n / [2\h sinh 2Xh 4- 4/( I 4- n>cosh 2Xh] sin A cos Ax </A/ A 

J 0 

/• -“50 

i xv (x.h) = 4p/n / [A/icoshA/i -(1 — n)/0 n) sinh A/ j] sinAsin\xdA/A 

Jo 


(65) 

(66) 

(67) 


Half-plane solutions 

The equations governing the behavior of the orthotropic half-plane are those for the 
corresponding strip, except that the conditions on the edge y = h are replaced by the re- 
quirement that the stresses vanish for x 2 + y 2 — ► sc in the upper half-plane. The Fourier 
sine and cosine transforms can again be usefully employed. 

Solution for type I material 

Because of the requirement that the stresses vanish for x 2 + y 2 -*■ oo, Eq. ( 1 8) can be 
replaced by 


<p c (A, y) = A(A)exp(-A*iy) 4- 5(A)exp(-Afc 2 v) (68) 


Satisfaction of the two transformed boundary conditions, Eqs. (12) and (13), leads to 


cp c (A, y) = /(A)/(fc, - k 2 )[k 2 exp(-Xkiy) - k i exp(-A* 2 y)] 

Then, the use of Eqs. (8M10) gives 

pOC 

a x = 2pk,k 2 /[n(k l - k 2 )l / [An exp(-A*,y) - k 2 exp(-A* 2 y)l sin A cos A.r d\/X 

Jo 

poo 

a v = 2p/[n(k\ - k 2 )] / [ki exp(-A k 2 y) - k 2 exp(-Afc,y)] sin A cos Ax d A/ A 

Jo 

pOC 

t, v = 2pk\k 2 /[n{ki - k 2 )] / [exp(-A* 2 y) - exp(-AAiy)] sin A sin Ax dA/A 

Jo 

All integrals can be evaluated in closed form [4]. The results are 

o T = 2 pk\k 2 /[n(k\ - Jfe 2 )](Ai{l/2 Arctan[2Aiy/(A?y 2 + * 2 ~ 01 + 5 i n / 2 ! 

— fc 2 {l/2Arctan[2fc 2 y/(A 2 y 2 4-x 2 - l)] + 5 2 7t/2}) (69) 

Oy = 2p/[n(ki — A 2 )](A| { 1/2 Arctan[2£ 2 y/(A 2 y 2 4-x 2 - 1)] 4- 5 2 7t/2} 

— £ 2 { 1/2 Arctan[2Aiy/(A 2 y 2 4- x 2 — 1)] 4- Si7t/2}) C^) 

T tv = pkik 2 /[2n(k\ — £ 2 )](ln{[A: 2 y~ + (x + l)"]/[^ 2 . v + (-* ~~ 0 ]} 

— ln{[fc 2 y 2 4- (x 4- O 2 ]/[0\v 2 4- (x - D 2 ]|) 


( 71 ) 


Transverse Stress Decay in a Specially Orthotropic Strip 


163 


where 


1 

k 2 y 2 

+ x 2 -l <0 

and 

S 2 = - 

h 

k\y 2 + x 2 — 

1 < 0 

0 

k 2 y 2 

+ x 2 - 1 > 0 J 


1° 

k\y 2 + x 2 - 

1 > 0 


and k[ and ki are the type I characteristic values defined earlier. 


(72) 


Solution for type III material 

For this material, after enforcement of the requirement that all stresses vanish for 
x 2 + y 2 ^ co in the upper half-plane, the transformed stress function takes the form 

cp c (A, y) = A(A)exp[-A(fci + ik 2 )y] + 5(A) exp[-A(*i — /A: 2 )v] (73) 

Satisfaction of the two transformed boundary conditions yields 

(p c (A, y) = — /(A)exp(— A/ciy)[cos \k 2 y + (^ 1 /^ 2 ) sin \k 2 y] (74) 


Then the use of Eqs. (8)— ( 1 0) gives 


/*0O 

c t x = —2p{k\ + k\)fn / [(k\/k 2 )smXk 2 y - cosAfc 2 y]exp(-A*iy)sinAcosA*dA/A 

Jo 

u y =2/7/71 1 [(/ci / k 2 ) sin \k 2 y + cos Xk 2 y] exp(— Afcj y) sin A cos A* <f A/A 

Jo 

P 00 

T xy — 2p(k 2 + k\)/{nk 2 ) / expt-M^sinA^sinAsinArdV^ 

Jo 

Again, all integrals can be evaluated in closed form [4], The results are 
a x = — p(k 2 + k^/n] i (k\/4k 2 )(^n\\k\y 2 + (1 + k 2 v + x) 2 \/\k]y 2 + (1 — k 2 y —■*)“]} 
+ ln{[^y 2 + (l + k 2 y - x) 2 ]/[k 2 { y 2 + (1 - k 2 y + x) 2 ]}) 

— 1/2 Arctan{2fciy/[fc 2 y 2 — 1 + ( k 2 y + x) 2 ]} — S\n/2 

— l/2Arctan{2fciy/[fc 2 y 2 — 1 + (k 2 y — x) 2 ]} — $2^/2] (75) 

<J y = p/n\(k\l4k 2 )(\n{\k\y 2 + (1 +k 2 y -\- x) 2 ]/[k\y 2 + (1 - k 2 y — x) 2 ]} 

+ ln{[fc 2 y 2 + (1 + k 2 y - x) 2 ]/[k]y 2 { 1 - k 2 y + x) 2 ]}) 

4- 1/2 Arctan{2&i,v/[fc 2 y 2 — 1 + (k 2 y + x) 2 ] J + Si7t/2 

+ 1/2 Arctan{2£|_v/[fc 2 }' 2 — 1 + ( k 2 y — x) 2 ]) + 527r/2] (76) 


where 


1 k\y 2 - 1 + (k 2 y + x) 2 < 0 

0 k 2 y 2 — 1 + (k 2 y + x) 2 > 0 


1 k\y 2 — 1 + (k 2 y - x) 2 < 0 1 
0 k\y 2 - 1 + (k 2 y - x) 1 > 0 j 

(77) 


and 

T, v = p(k 2 + k;) /(27ik 2 ){Arct<in[(k 2 y +x - l)/(*iy)] + Arctan[(* 2 :y -x+ l)/(*iy)] 
- ArctanK^y + x + l)/(*iy)] - Arctan[(^ 2 .v - x- l)/(*iy)]} (78) 
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Solution for isotropic material 

This solution is readily available in the literature [5]. The stresses are 

(j x = (p/7T){Arctan[ v/(.r - 1)] - Arctan[v/(.r + 1)] 

+ 2 y(jt 2 -r- i )/Kr + y 2 - l) 2 -h 4y 2 ]} (79) 

tJy = (p/7t){Arctan[y/{.r - 1)] - Arctan[v/(.r + 1)] 

- 2y(;t 2 - v 2 - l )/ [f.r 2 + y 2 - l) 2 + 4y 2 ] } (80) 

T,y = ( p/n)[4xy 2 /[(x 2 + y 2 - l) 2 +4y 2 ]} (81) 

where the Arctan functions have the range [0. 7t). 


Results and discussion 

In order to generate illustrative numerical and graphical results for stresses in some 
specially orthotropic configurations, a standard laminate analysis routine was used to 
compute equivalent orthotropic -sheet properties for a wide variety of symmetric angle- 
ply laminates in which each graphite-epoxy lamina had the following properties: E L = 
128 GPa, Et = 11.0 GPa, r| Lr =0.35, and Gu =5.74 GPa, where subscripts L and T 
denote longitudinal and transverse directions relative to the fibers. Table 1 contains a list of 
20 balanced and symmetric 24-ply laminates along with their equivalent orthotropic pro- 
perties, their corresponding characteristic values, and their designation as type I, type III, 
or the one quasi-isotropic type II. Ply angles are measured from the positive x axis. 


Table 1 

Elastic constants and characteristic values of various laminates (moduli given in GPa) 


Laminate 

£*(KT 6 ) 

£,(10- 6 ) 

O.tv 

GxvdO 6 ) 

Type 


k 2 

[0] 2 4 

128 

11.0 

0.350 

5.74 

I 

4.58 

0.742 

[±5] 6S 

125 

11.0 

0.421 

6.56 

I 

4.19 

0.803 

[±10] 6S 

118 

11.1 

0.617 

8.96 

I 

3.31 

0.983 

[±15]6s 

105 

11.2 

0.879 

12.5 

I 

2.12 

1.44 

[±20] 6S 

86.9 

11.4 

1.11 

17.0 

III 

1.45 

0.813 

[±25] 6S 

67.2 

11.8 

1.24 

21.8 

III 

1.16 

1.02 

[±30] 6S 

49.4 

12.5 

1.21 

26.3 

III 

0.924 

1.06 

l±35] 6S 

35.5 

13.9 

1.09 

29.9 

III 

0.745 

1.02 

[±40] 6S 

26.0 

16.1 

0.913 

32.3 

III 

0.616 

0.944 

[±45] 6S 

19.9 

19.9 

0.731 

33.1 

III 

0.534 

0.846 

[±50] 6S 

16.1 

26.0 

0.565 

32.3 

III 

0.485 

0.743 

[±55]6s 

13.9 

35.5 

0.425 

29.9 

III 

0.465 

0.639 

[±60]6s 

12.5 

49.4 

0.309 

26.3 

III 

0.467 

0.536 

[±65] 6S 

11.8 

67.2 

0.217 

21.8 

III 

0.487 

0.427 

[±70] 6S 

11.4 

86.9 

0.146 

17.0 

III 

0.525 

0.294 

[±75] 6S 

11.2 

105 

0.094 

12.5 

I 

0.472 

0.694 

[±80]6s 

1 1.1 

118 

0.058 

8.96 

I 

0.302 

1.02 

[±85] 6S 

11.0 

125 

0.037 

6.56 

I 

0.239 

1.25 

[90] 24 

11.0 

128 

0.030 

5.74 

I 

0.218 

1.35 

Isotropic 

50.9 

50.9 

0.312 

19.4 

II 

1 

1 
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In the numerical computation of stresses from either improper integrals or closed- 
form solutions, and in the graphical representation of the results, extensive use was made 
°f Mathematica™ software [6J. Since all of the improper integrals in problems A and B 
required numerical integration, considerable experimentation was performed to establish a 
practical upper limit in place of A = oo. An upper integration limit of A = 20 was found to 
be more than sufficient in all cases examined. 

The closed forms of the solutions to the orthotropic and isotropic half-plane problems 
make it convenient and, it is hoped, instructive to present some half-plane results in advance 
of results for the strips. 


Half-plane constant-stress contours 

To illustrate some of the effects that material properties can have on load diffusion 
rates and patterns, contour plots for the normalized stresses cr y /p and r xy /p are pre- 
sented in Figures 1 and 2, respectively, for the isotropic half-plane and for orthotropic 
half-planes which correspond to the following seven laminates: [OH, [±15] 6S , [±30] 6 s 
[±45] 6S , [±60] 6S , [±75] 6S , and [90] 24 . To aid in making direct comparisons, all plots have 
been generated for the same square region defined by (-6 < .r <6,0<y < 12}. 

Figure 1, which contains the five level curves u y /p = 0.9, 0.7, 0.5, 0.3, and 0.1 for 
each of the eight laminates, shows clearly that diffusion of the external loading into the 
half-plane is strongly influenced by laminate properties. For example, the maximum depth 
(i.e., distance, normalized by the width of the loaded region, from the edge y — 0) at which 
(Jy/p persists at the 10% level varies from about 6 for the [±30] 6S laminate to very much 
more than 12 for the [90] 2 4 laminate. The corresponding depth for the isotropic laminate 
lies well within these extremes. 

In addition, among the eight laminates there are large differences in the shapes of 
corresponding constant-cr v /p contours, especially for the lower levels of cr y /p. In fact, 
contours for the intermediate-ply-angle laminates, i.e., those composed of ±30°, ±45°, 
and ±60 c plies, exhibit a fundamental shape change as c j y /p becomes small. For large 
values of cr v / p the contours display an absolute maximum along the y axis, their axis of 
symmetry, whereas for small values of c J y /p they develop a relative minimum there. This 
behavior is consistent with the fact that the intermediate-angle-ply laminates have preferred 
directions of load propagation that are oblique to the direction of loading (90°), a feature 
that is especially evident in the [±45] 6S laminate (Figure le). 

Before contour plots of normalized shear stress, T xv /p, for the eight laminates were 
generated, ranges for r xy / p had to be established. This was done by examining the limiting 
behavior of T tv in the neighborhood of (.*, y) = (1, 0), a point of discontinuity of the 
edge loading, where the shear stress has an extreme value. For T max = Iim v ^o x* v (l, v), the 
following expressions were obtained: 


^ma x/P — 


*1*2 *l_ 

tc(£i — kf) ^ k\ 

k ? + k ; k^ 

— - Arc tan — 

nki k\ 

1 


71 


(type I) 
(type III) 
(isotropic) 



i 




< 


Upon appropriate passage to the limit, the type I and type III expressions reduce to the 
isotropic result. Numerical values of the maximum shear stresses in the eight featured 




In Figure 2, level curves are shown for t iv jp = ±0.05 and for equal increments of 
±0 05 to the I maximum } value, for which a level curve for x tv /p could be generated. Hence, 
the number of curves shown for a laminate is indicative of the seventy of the shear stress 
field in that laminate. Also, there is positive correlation between the severity of the shear 
stress field and the rate of decay of the direct stress due to the normal loading. For example, 
in the [±15]6s laminate the high shear stresses (Figure 2c) are consistent with the highly 
compact o y field (Figure lc). Similarly, in the [90] 24 laminate the benign shear stress field 
(Figure 2h) correlates with the gradual decay of a y (Figure \h). 

These stress contour plots for the half-plane have been presented for two reasons. First 
the relatively simple closed-form solutions for the half-plane lend themselves to economical 
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figure 2. Constant-x xy plots for eight different laminates. (Continued) 


ontour-plot generation, while the strip solutions, in their somewhat involved improper 
ntegral forms, do not. Second, although the stress contour plots for the half-plane cann 
v expected to mirror exactly the behavior of the elastic strips, the results may be helpfu 
u visualizing the corresponding responses of the strips, if proper account can be taken of 
.fleets of the support conditions along the edge y = h. Information from the half p 
Solutions that might be useful in this regard is y„, the 90% decay distance for a v , which 
:s defined here as the minimum (normalized) distance from the loaded edge for which 

j ix vh)/p < 0.1 for all x. 

' To find y 4 for most of the laminates, it sufficed to examine a v only along the y axis 
Vcause every level curve for ct v has a maximum there. However as was noted earlier, 
w the [±30] 6S and [±60] 6S laminates, and presumably for all balanced angle-ply lam 
nates between these, the level curves o,/p = 0.1 exhibit a relative minimum .on the 
, axis Thus, for those laminates the cjp = 0.1 level curves had to be examined fur- 
ivcr to find the correct value of v„. For the 20 laminates studied here, the values of y d 
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Figure 2. Constant- r, y plots for eight different laminates. (Continued). 


ranged from 6.05 to 33.7. These 90% decay distances for the half-planes are presented 
and discussed in more detail in connection with corresponding results for the elastic 

strips. 


Strip widths that ensure 90% decay 

A question that arises is whether the relatively convenient solutions for the half-plane 
can be useful in the design of similarly loaded strips. A partial answer might be obtained^ 
finding whether there is correlation between the 90% decay distances for the half-planes, 
y d and the minimum strip width required to ensure 90% decay of <r v across the strip. To 
this end, the formulas for given by Eqs. (24). (31) (37). <52 . (59. 

used to find h d , the minimum value of h for which a y (x , h)/p < 0. 1 for all x. The process 
which involved calculating tr v (x, h) in each laminate for numerous values of h. yielded 
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Table 2 

Comparisons of y d with minimum strip widths that 
ensure 90% decay of er v 



Laminate 

yj, Half-Plane 

[0124 

9.91 

[±5] 6S 

9.39 

[± 1 OJfts 

8.35 

[±15] 6S 

7.39 

[±20J 6S 

6.66 

[±25] 6S 

6.18 

[±30]6s 

6.05 

[±35] 6S 

6.51 

[±40] 6s 

7.30 

[±45] 6S 

8.25 

[±50] 6S 

9.27 

[±55] 6S 

10.4 

[±60] 6s 

12.0 

[±65] 6S 

14.8 

[±70] 6S 

18.4 

[±75] 6S 

22.6 

[±80] 6s 

27.2 

[±85] 6s 

31.5 

[90] 24 

33.7 

Isotropic 

12.7 


hj. Strip A h d , Strip B 


15.0 

14.6 

14.1 

13.6 

12.4 

11.6 

10.7 

9.54 

9.42 

7.74 

8.48 

6.37 

8.21 

6.49 

8.98 

7.43 

10.2 

8.70 

11.6 

9.96 

13.0 

11.1 

14.4 

11.9 

16.2 

12.8 

20.3 

15.2 

26.0 

21.4 

32.7 

29.1 

40.3 

37.8 

47.4 

45.7 

50.9 

49.5 

18.3 

16.2 



in Table 2 ’ Which also contains the va l»es of y d for the corresponding 

As was the case with some of the type III half-planes, the normal stress <r v (x h d ) on 
the supported edge of some type III strips has a relative minimum on its axis of symmetry 
whtch ,s the v axis. An obvious conclusion from Table 2 is that, for either strip 

ih , h UeS ° f kd ' 016 nUn,mUm Stnp Wldth re 9 uired to ens ure 90% decay of o v , are greater 
than the corresponding values of y d , the 90% decay distance in the half-plane. For the strip 

with the partially fixed edge (problem A), the ratio h d /y d varies only moderately amo^g 
the Inmates, from a low of 1.35 for the [±60 ]6S laminate to a high of 1.51 forLTo] 2 f 
and [90] 24 laminates. Thus, the use of a single h d ,y d ratio of abom 1.5 for 

reasonable 18 3 b orthotro P ,c stri P on die basis of the half-plane solution appears 

For the strip with the fixed edge (problem B), however, the range of h d /x d is consid- 

Xrt:- t a lo : of 1 03 for “* i±25i “ ■» ■ 4 <* lz ;»£ 

for i ] 24 nf h n ^s WI r St S ° f hd/yd falling beiow l m Thus, a single value 
d! d about 1.5 would again serve as an effective upper bound on the panel width 

required to ensure 90% decay of a,. However, because of the larger spread m C/Tthls 

approach could, for some laminates, lead to wider than necessary strips. Also, in the case 

amrie Tab ' e 2 ,nd,cates a stro "g interaction between the intermediate-ply 

angles and the restraint against tangential displacement along the supported edge, which 
can result in substantial shear stresses there. 
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°y(K.h d > 06 
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{ Isotropic, h d = 18.3 
C± 1S 1 6 .. h,, = 10.7 
[± 7 5W h d = 32.7 


-/f°l 2 4. ^*15.0 
\m 24 , h d = 50.9 


<yx,h d ) 06 

P .04 


±30]*, h d = 8.21 
± ®0] 6s . h d = 16.2 


-[±45W h d *11.6 



a) Isotropic and four type-1 laminates 


b) Three type-1 II laminates. 


Figure 3. Graphs ofo y (x. h d )/p for eight laminates (problem A). 

Normal stress on the supported edge 

To illustrate the variety of normal stress distributions that occur along the supported 
ge in response to the localized uniform loading on the opposite edge, af.x.h d )/p is 
p otted as a function ot ,t in Figure 3 for problem A (partially fixed edge) and in Figure 4 

thar P H° b , e , m h edge) ’ f ° r thC S3me dght ,aminates featured e^lier. It is worth noting 

fnH qm y r-c°si 7? graphS ° f CTvU ’ hd)/p for the Allowing laminate pairs: [OK 
and [90] 24 [± 5] 6S and [±75]«. and [±30] 6S and [±60] 6S . These identities result from 
a reciprocal relationship that exists between characteristic values of these pair members 
because one pair member is the 90' rotation of the other. Note, however, that although the 
SscSiff SUPPOrtCd Cdge ° f these reci Procal pairs are identical (except forpossible 

S l d t renCCS ? t0 nUmenCal r ° Undoff) ’ they a PP ly to ver y different strip widths, 
cidentally the reason for the rather remarkable near-identity in both problems between 

^ graphs of <Jy(x. hj)/p for the isotropic material and for the reciprocal pair [±15kc and 
L± /5]6 s is not apparent to the author. Jfti> 

F | g !! reS 3 l and Sh ° W th3t graphS ° f CTv(r ’ hd)/p for the i^tropic and type I materials 
! ‘ either qualitatively or quantitatively, between problem A and problem B. Figures 
3h and 4 b show that graphs of a v (x, h d )/p for the type III materials are also qualitatively 
similar m problems A and B; however, their relative minima on their line of symmetry are 



( Isotropic, h d = 16.2 
[-^ 5]^, h d = 9.54 
[±75]^, h d = 29.1 

. 4 [0J24. h d = 14 

V l[90] 2 „h d = 4! 


<T y (x,h d ) ' 
p .04 


a) Isotropic and tour type-1 laminates 


AK.: K _ rvr- 


301s,, h d = 6.49 
Wle*. h d = 12.8 



Figure 4. Graphs oftj y (x, hi) /Pfor eight laminates ( problem B). 


b) Three type-ill laminates. 
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u-v=0 

( 

( J 



Figure 5. Graphs of r xy (x, h^)/p for five type UI laminates. 


more pronounced in problem B, due primarily to the tangential restraint along the edge 
v = h d . In any event, if a design requirement is at least 90% decay of o\ across the strip, 
then for a large number of angle-ply laminates. Table 2 identifies the minimum width needed 
to satisfy that requirement. 

In the case of problem B, however, there is the possibility of significant shear stresses 
on the fixed edge y = h d . With that possibility in mind, shear stresses along y = h d were 
calculated using Eqs. (53), (60), and (67) for type I, type III. and isotropic material, respec- 
tively. Significant shear stresses (defined here as satisfying T xy /p > 0. 1) on the edge y = h d 
were not found in any of the isotropic or type I materials studied. In the type III materials, 
however, shear stresses that meet the significance criterion were present in the following 
five laminates: [±25] 6S , [±30] 6S . [±35] 6S , [±40] 6S , and [±45] 6S . Plots of T, y (jc f h d )/p for 
these five laminates are shown in Figure 5. By a slight margin, maximum shear stresses on 
the supported edge are greatest in the [±30]6s and [±35]6s laminates. If shear stresses of 
this magnitude are of little concern, then the list of h d values given in Table 2 should suffice 
as minimum panel widths for design purposes. On the other hand, if these shear stresses are 
unacceptably large, then Eq. (60) can be used along with appropriate numerical integration 
tools to determine new, greater panel widths that ensure sufficient decay of both cr v and T tv . 
In the case of the [±30]6$ laminate, for example, this augmented criterion would require 
that the minimum strip width be increased from 6.49 to 7.61. 

In the case of problem B, an additional potential complication exists. As a result of 
the fixed boundary condition, the relation a x =T} xv <j y holds everywhere along the edge 
v = /j, so that in the small group of laminates for which rj t v > 1 (see Table 1), it follows 
that <j t > <t v there. Thus, if the 90% decay criterion were to be applied to c x as well as 
to cr v , then for r[ xy > 1 the condition o\(.t, h d )/p < 0.1 would have to be replaced by the 
condition ri xv cr v (jt, h d )/p <0.1, which would result in greater minimum strip widths for a 
few laminates. 

A trend that might, upon first consideration, appear to be counterintuitive can be seen in 
Table 2. Note that for every laminate the strip width that ensures 90% decay of a v (.t, h d )/p 
is greater in problem A (partially clamped edge) than in problem B (clamped edge). This 
condition is equivalent to the following: for identical laminates of equal strip width, the 
normal boundary stress that reacts the applied edge loading, i.e., o\(.y, h ), has a greater 
maximum value on the partially fixed edge than on the "stiffer*’ fully fixed edge. Such a trend 
runs counter to the popular adage that suffer boundaries, by offering greater resistance, incur 
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a) Problem A - partially fixed edge. 


b) Problem B - fully fixed edge. 


Figure 6. Shear stress contour plots for the isotropic cases of problems A and B with h = 5. 


greater stress. This apparent anomaly can be resolved by referring to Figures 6 a and 6 b y 
which contain shear stress contour plots for the isotropic cases of problems A and B, 
respectively, with normalized strip width arbitrarily taken as 5. (The same point can also 
be illustrated with any of the other laminates and any small-to-moderate strip width.) Note 
the considerably greater extent of the shear stress field in problem B, as well as the nonzero 
shear stress distribution on the fixed edge in problem B. At all points on the partially fixed 
edge of problem A, of course, shear stress is required to be zero. In effect, the restraint 
against tangential displacement along the fixed edge gives rise to the more extensive shear 
stress field by which the broader diffusion of the applied normal load is accomplished. 
A consequence of this broader diffusion of cr v is lower maximum normal stress on the 
supported edge. In other words, by promoting greater diffusion of cr v , the fixed boundary 
can be said to react the applied load more efficiently than does the partially fixed boundary. 

Concluding remarks 

Solutions in the form of improper integrals have been presented for the stresses in a 
specially orthotropic infinite strip with one edge partially fixed or completely fixed while 
the other edge is subjected to localized uniform normal loading. In addition, closed-form 
solutions for the stresses in a similarly loaded specially orthotropic half-plane are obtained. 
The results are used to generate illustrative contour plots of the direct stress and the shear 
stress in various half-planes, and to determine minimum strip widths required to ensure 
90% decay of stress across the strip. 

While the graphical and tabular results presented are for specific laminated composite 
materials, the formulas for the solutions can be used to calculate the stresses in any similarly 
loaded specially orthotropic laminated strip or half-plane whenever its elastic constants are 
known. 
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Appendix: Solutions for isotropic strip problems 

The method of solution for the isotropic strip problems is identical to that for the 
orthotropic strip problems; however, in the isotropic case the algebra is somewhat simpler. 
For isotropic material, Eqs. (1) become 






y.xy = “ y +n,jc 


I 

G 


(Al) 


Equations (2) still give the stresses in terms of the Airy stress function; however, Eq. (3) 
simplifies to 


+ 2(p XX yy + tp WW — 0 (A2) 

Fourier transforms [see Eqs. (8) and (9)] are again employed. Equation (10) becomes 
= (PC.yy C = “A 2 (pc TtvS = Acpc.v 

(A3) 

U S = (\/\E)(ipc,yy -F TjiV’Cpc) Vc = (l/X 2 E)((p c ,yyy ~ (2 + r|)A 2 (pc, v ) 
and transformation of Eq. (A2) gives 

VC.yy.vy — 2A~(pc. vv “F A 4 (pc = 0 (A4) 

For both strip problems, the transformed stress function has the form 

<Pc(A,y) = [A(A) + y£(A)]coshAy + [C(A) + yD(A)] sinh Ay (A5) 


Problem A: Partially fixed upper edge 

The transformed boundary conditions are given by Eqs. ( 1 2)— ( 15). Substitution of 
Eq. (A5) into the transformed boundary conditions leads to four simultaneous equations 
which have the solution 


A = -/(A) 

B = — 2A 2 /(A) sinh 2 A/z/ A 
C — 2A /(A) sinh 2 A/t / A 
D = 2A 2 /(A)sinh \h cosh A/z/ A 


where 


A = A(sinh 2A/z + 2A h) 
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Substitution of Eq. (A6) into Eq. (A5), followed by use of the first three of Eq. (10) and 
Fourier inversion, leads to the stress formulas given by Eqs. (32M37). 

Problem B: Fixed upper edge 

The transformed boundary conditions are given by Eqs. (52)— <55). Substitu- 
tion of Eq. (A5) into Eqs. (52)— (55) leads to simultaneous equations which have the 
solution 

A = -/(A) 

B = -A/(A)[(3 - r|) sinh \h cosh \h - (1 + q)A/i]/D 

C = A/(A)[(3 — q) sinh A/i cosh Xh - ( 1 -f- r\)\h]/D 

D = A/(A)[(3 — q) sinh 2 Xh + 2]/D (A7) 

where 


D = /[(3 - n) sinh 2 Ih + (1 + n)l 2 h 2 + 4/(1 +«)] 

Substitution of Eq. (A7) into Eq. (A5), use of the first three of Eq. ( 10) and Fourier inversion 
yield the stress formulas given by Eqs. (61M67). 


